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ABSTRACT 


The interaction of a runaway supermassive black hole (SMBH) with the circumgalactic medium (CGM) can 
lead to the formation of a wake of shocked gas and young stars behind it. Here we report the serendipitous 
discovery of an extremely narrow linear feature in HST/ACS images that may be an example of such a wake. 
The feature extends 62 kpc from the nucleus of a compact star-forming galaxy at z = 0.964. Keck LRIS spectra 
show that the [O ri1]/H ratio varies from ~ 1 to ~ 10 along the feature, indicating a mixture of star formation 
and fast shocks. The feature terminates in a bright [O 111] knot with a luminosity of zz 1.9 x 10*! ergs s™!. The 
stellar continuum colors vary along the feature, and are well-fit by a simple model that has a monotonically 
increasing age with distance from the tip. The line ratios, colors, and the overall morphology are consistent with 
an ejected SMBH moving through the CGM at high speed while triggering star formation. The best-fit time 
since ejection is ~ 39 Myr and the implied velocity is vgy ~ 1600 kms !. The feature is not perfectly straight 
in the HST images, and we show that the amplitude of the observed spatial variations is consistent with the 
runaway SMBH interpretation. Opposite the primary wake is a fainter and shorter feature, marginally detected 
in [O nI] and the rest-frame far-ultraviolet. This feature may be shocked gas behind a binary SMBH that was 


ejected at the same time as the SMBH that produced the primary wake. 


1. INTRODUCTION 


There are several ways for a supermassive black hole 
(SMBH) to escape from the center of a galaxy. The first step 
is always a galaxy merger, which leads to the formation of 
a binary SMBH at the center of the merger remnant (Begel- 
man et al. 1980; Milosavljevié & Merritt 2001). The binary 
can be long-lived, of order ~ 10? yr, and if a third SMBH 
reaches the center of the galaxy before the binary merges, a 
three-body interaction can impart a large velocity to one of 
the SMBHs leading to its escape from the nucleus (Saslaw 
et al. 1974; Volonteri et al. 2003; Hoffman & Loeb 2007). 
Even in the absence of a third SMBH, the eventual merger 
of the binary can impart a kick to the newly formed black 
hole through gravitational radiation recoil (Bekenstein 1973; 
Campanelli et al. 2007). The velocity of the ejected SMBH 
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depends on the mechanism and the specific dynamics. Gen- 
erally the kicks are expected to be higher for slingshot scenar- 
ios than for recoils (see, e.g., Hoffman & Loeb 2007; Kesden 
et al. 2010), although in exceptional cases recoils may reach 
~ 5000 km s^! (Campanelli et al. 2007; Lousto & Zlochower 
2011). In both scenarios the velocity of the SMBH may ex- 
ceed the escape velocity of the host galaxy (see, e.g., Saslaw 
et al. 1974; Hoffman & Loeb 2007; Lousto et al. 2012; Ri- 
carte et al. 2021b). 

Identifying such runaway SMBHs is of obvious interest but 
difficult. The main focus has been on the special case where 
the black hole is accreting at a high enough rate to be iden- 
tified as a kinematically or spatially displaced active galac- 
tic nucleus (AGN) (Bonning et al. 2007; Blecha et al. 2011; 
Komossa 2012). For such objects, the presence of a SMBH 
is not in doubt, but it can be difficult to determine whether 
they are “naked” black holes or the nuclei of merging galax- 
ies (see, e.g., Merritt et al. 2006). Candidates include the 
peculiar double X-ray source CID-42 in the COSMOS field 


2 


(Civano et al. 2010) and the quasars HE0450-2958 (Mag- 
ain et al. 2005), SDSSJ092742943 (Komossa et al. 2008), 
E1821+643 (Robinson et al. 2010; Jadhav et al. 2021), and 
3C 186 (Chiaberge et al. 2017). 

Quiescent (non-accreting) runaway SMBHs can be de- 
tected through the effect they have on their surroundings. 
As noted by Boylan-Kolchin et al. (2004) and discussed in- 
depth by Merritt et al. (2009), some of the stars in the nuclear 
regions of the galaxy are expected to remain bound to the 
SMBH during and after its departure. The stellar mass that 
accompanies the black hole is a steeply declining function 
of its velocity, and generally < Mgg. This leads to peculiar 
objects, dubbed "hyper compact stellar systems" (HCSS) by 
Merritt et al. (2009), with the sizes and luminosities of glob- 
ular clusters or ultra compact dwarf galaxies but the velocity 
dispersions of massive galaxy nuclei. HCSSs could there- 
fore be easily identified by their kinematics, but measuring 
velocity dispersions of such faint objects is difficult beyond 
the very local Universe. Other potential detection methods 
include gravitational lensing (Sahu et al. 2022) and tidal dis- 
ruption events (e.g., Ricarte et al. 2021a; Angus et al. 2022). 
No convincing candidates have been found so far. 

Another way to identify runaway SMBHs is through the 
effect of their passage on the surrounding gas. This topic 
has an interesting history as it is rooted in AGN models that 
turned out to be dead ends. Saslaw & De Young (1972) in- 
vestigated the suggestion by Burbidge et al. (1971) and Arp 
(1972) that the redshifts of quasars are not cosmological but 
that they were ejected from nearby galaxies. In that context 
they studied what happens when a SMBH travels supersoni- 
cally through ionized hydrogen, finding that this produces a 
shock front with a long wake behind it. Shocked gas clouds 
in the wake can cool and form stars, potentially illuminating 
the wake with ionizing radiation from O stars. Rees & Saslaw 
(1975) analyzed the possibility that double radio sources are 
produced by the interaction of escaped SMBHs with the in- 
tergalactic gas. They find that this is plausible from an ener- 
getics standpoint, although now we know that the alternative 
model, feeding of the lobes by jets emanating from the nu- 
cleus (Blandford & Rees 1974), 1s the correct one. 

Perhaps because of these somewhat inauspicious connec- 
tions with failed AGN models there has not been a great deal 
of follow-up work in this area. To our knowledge, the only 
study of the formation of wakes behind runaway SMBHs in a 
modern context is de la Fuente Marcos & de la Fuente Mar- 
cos (2008), who analyze the gravitational effect of the pas- 
sage of a SMBH using the impulse approximation. They find 
that the SMBH can impart a velocity of a few to several tens 
of km s^! on nearby gas clouds, and that the gas can then be- 
come unstable to fragmentation and star formation. The out- 
come is qualitatively similar to the analysis of Saslaw & De 


Young (1972), in the sense that, under the right conditions, 
star formation can occur along the path of the SMBH. 

In this paper we report on the serendipitous discovery of 
a remarkable linear feature in HST images that we suggest 
may represent such a SMBH-induced wake. We also identify 
two candidate hyper-compact stellar systems, one embedded 
in the tip of the wake and the other on the opposite side of 
the galaxy from which they may have escaped. 


2. A 62 KPC LONG LINEAR FEATURE AT Z = 0.964 
2.1. Identification in HST/ACS images 


We serendipitously identified a thin, linear feature in HST 
ACS images of the nearby dwarf galaxy RCP28 (Román 
et al. 2021; van Dokkum et al. 2022a), as shown in Fig. 
1. RCP28 was observed September 5 2022 for one orbit in 
F606W and one orbit in F814W, in the context of mid-cycle 
program GO-16912. The individual £1c files were com- 
bined using DrizzlePac after applying a flat field correc- 
tion to account for drifts in the sensitivity of the ACS CCDs 
(see van Dokkum et al. 2022b). Upon reducing the data an 
almost-straight thin streak was readily apparent in a visual as- 
sessment of the data quality (see Fig. 1). Based on its appear- 
ance we initially thought that it was a poorly-removed cos- 
mic ray, but the presence of the feature in both filters quickly 
ruled out that explanation. The total AB magnitude of the 
streak is F814W = 22.87 + 0.10 and its luminosity-weighted 
mean color is F606W —F814W = 0.83 + 0.05. 

The streak points to the center of a somewhat irregular- 
looking galaxy, at a = 2541™45843; 6 = —8?20'55"4 (32000). 
The galaxy has F814W = 21.86 + 0.10 and F606W — 
F814W = 0.84 + 0.05; that is, the brightness of the streak is 
& 40 96 of the brightness of the galaxy and both objects have 
the same color within the errors. Not having encountered 
something quite like this before in our own images or in the 
literature, we decided to include the feature in the observing 
plan for a scheduled Keck run. 


2.2. Redshift 


The feature was observed with the Low-Resolution Imag- 
ing Spectrometer (LRIS; Oke et al. 1995) on the Keck I tele- 
scope on October 1 2022. The 300 lines mm ! grism blazed 
at 5000 A was used on the blue side and the 400 lines mm^! 
grating blazed at 8500 À on the red side, with the 680nm 
dichroic. The 170 longslit was used, centered on the galaxy 
coordinates with a position angle of 327°. The total exposure 
time was 1800s, split in two exposures of 900 s. Conditions 
were good and the seeing was ~ 08. On October 3 we ob- 
tained a high resolution spectrum with the 1200 lines mm! 
grating blazed at 9000 À in the red. Five exposures were ob- 
tained for a total exposure time of 2665 s. Conditions were 
highly variable, with fog and clouds hampering the observa- 
tions. 
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Figure 1. Top left: F606W + F814W HST/ACS image of the linear feature and its surroundings. Top right: Zoomed view of the F606W image. 
The feature shows a compact bright spot at the narrow tip, and seems to broaden toward the galaxy. Bottom left: Color image, generated from 
the F606W and F814W images. Bottom right panels: Sections of LRIS spectra near bright emission lines. The feature and the galaxy are at the 
same redshift. The kinematics and line strengths show complex variations along the feature. 


Data reduction followed standard procedures for long slit 
observations. Sky subtraction and initial wavelength calibra- 
tion were done with the PypeIt package (Prochaska et al. 
2020). The wavelength calibration was tweaked using sky 
emission lines, and the data from the individual exposures 
were combined. A noise model was created and cosmic rays 
were identified as extreme positive deviations from the ex- 
pected noise. For the low resolution spectrum a relative flux 
calibration, enabling the measurement of line ratios, was per- 
formed using the spectrophotometric standard HS 2027. 

We find continuum and strong emission lines associ- 
ated with the feature. The lines are readily identified as 
the redshifted [O 11] A43726,3729 doublet, Hy, H8, and 
[O 111] A44959,5007. The redshift is z = 0.964, and the im- 
plied physical extent of the feature, from the nucleus of the 
galaxy to its tip, is 62 kpc. The 2D spectrum in the regions 
around the strongest emission lines is shown in the bottom 
panels of Fig. 1. The lines can be traced along the entire 
length of the feature. There are strong variations in the line 
strengths and line ratios, as well as in the line-of-sight ve- 
locity. We will return to this in following sections. The S/N 
ratio in the high resolution spectrum is low, about 1/4 of that 
in the low resolution spectrum. 


3. PROPERTIES OF THE HOST GALAXY 
3.1. Morphology 


The same emission lines are detected in the galaxy, con- 
firming that it is at the same redshift as the linear feature (see 
Fig. 1). The galaxy is compact and somewhat irregular, as 
shown in Fig. 1 and by the contours in Fig. 2. We determine 
the half-light radius of the galaxy with galfit (Peng et al. 
2002), fitting a 2D Sersic profile and using a star in the image 
to model the point spread function. We find re ~ 1.2 kpc, but 
we caution that the fit has significant residuals. The irregu- 
lar morphology may be due to a recent merger or accretion 
event, although deeper data are needed to confirm this. 


3.2. Ionization mechanism 


We measure the strength of the strongest emission lines 
from the 2D spectra. The continuum was subtracted by fitting 
a first-order polynomial in the wavelength direction at all spa- 
tial positions, masking the lines and their immediate vicin- 
ity. Line fluxes were measured by doing aperture photome- 
try on the residual spectra. No corrections for slit losses or 
underlying absorption are applied. We find an [O 111] flux of 
F-(10ED x10 ergs s! cm? and [O IH = 1.90.2. 

The interpretation of the line fluxes depends on the ioniza- 
tion mechanism, which can be determined from the combina- 
tion of [O rrr]/H8 and [N rr]/Ho. Ha and [N 1] are redshifted 
into the J band, and we observed the galaxy with the Near- 
Infrared Echellette Spectrometer (NIRES) on Keck II on Oc- 
tober 4 2022 to measure these lines. NIRES provides cross- 
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Figure 2. Morphology of the galaxy in F606W and F814W. The 
arrow indicates the direction of the linear feature. The galaxy is 
compact, with a half-light radius of re = 1.2 kpc, and shows irregular 
features possibly indicating a recent merger and/or a connection to 
the linear feature. 


dispersed near-IR spectra from 0.94 um — 2.45 um through a 
fixed 0755 x 18" slit. A single 450 s exposure was obtained 
in good conditions, as well as two adjacent empty field expo- 
sures. In the data reduction, the empty field exposures were 
used for sky subtraction and sky lines were used for wave- 
length calibration. The Ha and [N r1] 46583 emission lines 
of the galaxy are clearly detected, as shown in the inset of 
Fig. 3. The emission lines of the galaxy are modeled with 
the redshift, the Ho line strength, the [N IT] line strength, and 
the velocity dispersion as free parameters. The best-fitting 
model is shown in red in Fig. 3. We find a velocity disper- 
sion of aga =60+7 km s^! and [N rr]/Ha = 0.23 +0.06, with 
the uncertainties determined from bootstrapping. The im- 
plied metallicity, using the Curti et al. (2017) calibration, is 
Z «-0.08*005. 

The location of the galaxy in the BPT diagram (Baldwin 
et al. 1981) is shown in Fig. 3. For reference, data from the 
Sloan Digital Sky Survey (SDSS) DR7 are shown in grey 
(Brinchmann et al. 2004). The galaxy is slightly offset from 
the SDSS relation of star-forming galaxies and quite far from 
the AGN region in the upper right of the diagram. The off- 
set is consistent with the known changes in the ISM condi- 
tions of star forming galaxies with redshift (see, e.g., Steidel 
et al. 2014; Shapley et al. 2015). The lines in Fig. 3 show 
the redshift-dependent Kewley et al. (2013) division beyond 
which AGN begin to contribute to the line ratios. The galaxy 
is well within the “pure” star formation region for z= 1. 


3.3. Star formation rate and stellar mass 


We infer the star formation rate of the galaxy from the HG 
luminosity, which is Lyg = (2.5 +0.5) x 10“! ergss !. The 
Kennicutt (1998) relation implies an approximate star forma- 
tion rate of 6 Mc yr ! for the dust-free case and 14 Mc yr! 
for 1 mag of extinction. The stellar mass of the galaxy can 
be estimated from its luminosity and color. We generate 
predicted F606W — F814W colors for stellar populations at 
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Figure 3. The location of the galaxy in the BPT diagram, with 
SDSS galaxies in light grey. The lines divide “pure” star forming 
galaxies from those with an AGN contribution to their line ratios, 
for z = 0 and z= 1 (Kewley et al. 2013). The location is as expected 
for az=1 star forming galaxy. The inset shows the NIRES spectrum 
in the Ha region. The red line is the best fit. 


z= 0.964 with the Python-FSPS stellar population model- 
ing suite (Conroy et al. 2009). We find that the observed color 
of the galaxy can be reproduced with a luminosity-weighted 
age of ~ 150 Myr and no dust or an age of ~ 65 Myr with 
Ay ~ 1. The implied stellar mass is Mga ~ 7 x 10° Mo. The 
typical star formation rate of a galaxy of this mass at z = 1 is 
= 8Mo yr ! (Whitaker et al. 2014), similar to the observed 
star formation rate. 

We conclude that the galaxy has normal line ratios and 
a normal specific star formation rate for its redshift. Its 
age is highly uncertain given that the color is dominated 
by the most recent star formation, but if we take the ~ 
100 Myr at face value, the past-average star formation rate is 
c- 70 Mc yr !, an order of magnitude larger than the current 
value. The galaxy shows morphological irregularities and is 
overall quite compact. Its half-light radius of 1.2 kpc is a fac- 
tor of ~ 3 smaller than typical galaxies of its stellar mass 
and redshift (van der Wel et al. 2014), which implies that its 
star formation rate surface density is an order of magnitude 
higher. Taken together, these results suggest that the galaxy 
experienced a recent merger or accretion event that led to the 
funneling of gas into the center and a burst of star formation 
~ 10? yr ago. 


4. SHOCKS AND STAR FORMATION ALONG THE 
FEATURE 


4.1. Variation in continuum emission and line ratios 


The linear feature is not uniform in either continuum 
brightness, color, line strengths, or line ratios. The variation 
along the feature in the F606W (Arest = 0.31 um) continuum, 
the F606W — F814W color, and in the [O 111] and H£ lines is 
shown in Fig. 4. Note that the spatial resolution of the contin- 
uum emission is ~ 8x higher than that of the line emission. 
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Figure 4. The four panels correspond to the rest-frame near- 
UV continuum, F606W — F814W color, [O ri], and H8 emission 
along the linear feature (pictured at the top). The F606W contin- 
uum shows strong variation on all spatial scales, and is brightest 
at the furthest point from the galaxy. The color shows large and 
seemingly random variations. The [O r11]/H£ ratio varies by a fac- 
tor of ~ 10 along the feature, with some regions likely dominated 
by shock ionization and others dominated by H II regions. 


There is a general trend of the continuum emission becom- 
ing brighter with increasing distance from the galaxy. The 
continuum reaches its peak in a compact knot at the tip; be- 
yond that point the emission abruptly stops. As shown in 
Fig. 1 the continuum knot at the tip coincides with a lumi- 
nous [O 11] knot in the spectrum. The [O rir] A5007 flux of 
the knot is F ~ 3.9 x 10^ ergs s™! cm™, and the luminos- 
ity is L ~ 1.9 x 10*! ergss^!. The [O rr1]/H ratio reaches 
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~ 10 just behind the knot, higher than can be explained by 
photoionization in H II regions. 

The ionization source could be an AGN, although as dis- 
cussed in more detail in § 6.4.3 the [OTII] emission is so 
bright that an accompanying X-ray detection might be ex- 
pected in existing Chandra data. An alternative interpreta- 
tion is that the bright [O 11] knot is caused by a strong shock 
(see Shull & McKee 1979; Dopita & Sutherland 1995; Allen 
et al. 2008). In the models of Allen et al. (2008), photoion- 
ization ahead of a fast (> 500km s7!) shock is capable of 
producing [O ri1]/Hf ~ 10, and the expected associated soft 
X-ray emission (Dopita & Sutherland 1996; Wilson & Ray- 
mond 1999) may be below current detection limits. There is 
at least one more region with elevated [O r11]/H5 ratios (at 
r ~ 25 kpc), and the [O III] emission near the tip could sim- 
ply be the strongest of a series of fast shocks along the length 
of the feature. 


4.2. Stellar populations 


In between the two main shocks is a region where O stars 
are probably the dominant source of ionization. At distances 
of 40 kpc < r < 50 kpc from the galaxy the [O r11]/Hf ratio 
is in the 1 —2 range and there are several bright continuum 
knots. These knots show strong F606W —F814W color varia- 
tion, mirroring the striking overall variation along the feature 
that was seen in Fig. 4. In Fig. 5 we compare the measured 
colors of three knots to predictions of stellar population syn- 
thesis models. They were chosen because they span most 
of the observed color range along the feature. The models 
span a metallicity range of —1 < Z < 0 and have either no 
dust (blue) or Ay 2 1 mag (red). The metallicity range en- 
compasses that of the galaxy (Z ~ —0.1). 

We find that the knots can indeed be young enough (< 
10 Myr) to produce ionizing radiation. However, it is difficult 
to derive any quantitative constraints as there is no straight- 
forward relation between age and color in this regime. The 
reason for the complex model behavior in Fig. 5 is that the 
ratio of red to blue supergiants changes rapidly at very young 
ages (“blue loops"; see, e.g., Walmswell et al. 2015). We note 
that the evolution of supergiants is uncertain (see, e.g., Chun 
et al. 2018) and while the overall trends in the models are 
likely correct, the detailed behavior at specific ages should 
be interpreted with caution (see, e.g., Levesque et al. 2005; 
Choi et al. 2016; Eldridge et al. 2017). In $ 7.1 we interpret 
the overall trend of the color with position along the feature 
in the context of our proposed model for the entire system. 

Finally, we note that the knots appear to have a character- 
istic separation, as can be seen in Fig. 5 and in the pattern of 
peaks and valleys from r = 30 kpc to r = 50 kpc in the F606W 
emission in Fig. 4. The separation is zz 4 kpc. This could be 
coincidence or be an imprint of a periodicity in the cooling 
cascade of the shocks. 
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Figure 5. Comparison of the observed colors of several knots in 
the feature (shown at the top) to model predictions of Conroy et al. 
(2009) for different ages. Dashed model predictions are for a metal- 
licty Z = —1, solid for Z = —0.5, and dot-dashed lines are for Z = 0. 
Blue lines are dust-free models and red lines illustrate the effect of 
dust attenuation with Ay = 1. Horizontal lines are measurements for 
the three knots. The ages of the youngest stars are likely < 30 Myr, 
but there is no straightforward relation between age and color in this 
regime. The observed colors span a similar range as the models and 
are consistent with a wide range of possible metallicities, ages, and 
dust content. 


5. A “COUNTER” LINEAR FEATURE ON THE OTHER 
SIDE OF THE GALAXY 


The LRIS slit covered the galaxy and the feature and also 
extended beyond the galaxy on the other side. There is no 
spatially-extended F606W or F814W emission on this side 
but there is an unresolved object, “B”, that is located at a 
distance of 4"4 from the galaxy within a few degrees of the 
orientation of the feature (see Fig. 6). The LRIS spectrum 
in the vicinity of the redshifted [O III] line is shown in the 
middle panel of Fig. 6, after subtracting the continuum and 
dividing by a noise model to reduce the visual effect of sky 
residuals. 

We detect a knot of [O11] A5007 emission near the lo- 
cation of B, redshifted by ~ 40kms"! with respect to the 
galaxy. Furthermore, there is evidence for faint [O 111] emis- 
sion in between the galaxy and B. This “counter” linear fea- 
ture is also seen in a u-band image, shown in the right panel 
of Fig. 6. The object was serendipitously observed with 
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Figure 6. Left: Section of the summed ACS F606W+F814W image, with the LRIS slit indicated in blue. Besides the tip of the linear 
feature, A, there are two other bright spots in the vicinity, B and C. Object B falls in the slit. Center: Section of the LRIS spectrum around the 
[O 111] A5007 line. Object B is detected, as well as faint emission in between B and the galaxy. The attached panel shows the intensity along the 
feature, on a logarithmic scale. Right: The presence of a “counter” feature is confirmed through its detection in the u-band, which samples the 
rest-frame far-UV. For clarity the u-band image was binned by a factor of 6 in the direction perpendicular to the slit (and then expanded again 
to retain the correct spatial scale). Also note that the primary feature extends all the way to the galaxy, in marked contrast to the pronounced 


gap between the galaxy and the feature in the ACS image. 


MegaPrime on the Canada France Hawaii Telescope (CFHT) 
on September 11 and 12, 2020 in the context of program 
20BO44 (PI: A. Ferguson). The total exposure time was 
11,880 s; the data reduction is described in M. L. Buzzo et 
al., in preparation. 

The u-band surface brightness of the counter feature is 
approximately 5x fainter than on the other side, and it ap- 
pears to terminate at the location of the [O 111] knot. Further- 
more, the primary feature extends all the way to the galaxy 
in the u-band: there is no gap at r < 25 kpc as is the case 
in the ACS data. The u-band samples the rest-frame far-UV 
(Arest % 0.18 um), and we conclude that the far-UV emission 
of the entire system is largely decoupled from the near-UV 
emission that is sampled with ACS. The total far-UV bright- 
ness of the linear emission is ~ 70 96 of the far-UV bright- 
ness of the galaxy, whereas this fraction is only ~ 40% at 
Arest © 0.36 um. 

The detection of the counter feature in the rest-frame far- 
UV shows that the [O 111] emission is likely real and caused 
by shocks. The combination of [O 111] line emission and far- 
UV continuum emission has been linked to cooling radiation 
of fast (= 100km s^!) shocks, both theoretically (e.g., Suther- 
land et al. 1993), and observationally, for instance in sections 
of supernova remnants (Fesen et al. 2021). 

It is difficult to determine the relationship between object 
B and the counter feature. It has F814W = 25.28 +0.10 (AB) 
and F606W —F814W = 0.84 + 0.14, and it is misaligned by 


4? from the line through A and the galaxy. We will discuss 
the nature of B in the context of our preferred overall model 
for the system in § 6.4. There is also another compact object, 
C, that is nearly exactly opposite to B in angle and distance. 
This object was not covered by the LRIS slit and we have no 
information about it, except that it is bluer than B. 


6. INTERPRETATION 
6.1. Various straight-line extragalactic objects 


With the basic observational results in hand we can con- 
sider possible explanations. Thin, straight optical features 
that extend over several tens of kpc have been seen before in 
a variety of contexts. These include straight arcs, such as the 
one in Abell 2390 (Pello et al. 1991); one-sided tidal tails, 
with the Tadpole galaxy (Arp 188) being the prototype (Tran 
et al. 2003); debris from disrupted dwarf galaxies, like the 
multiple linear features associated with NGC 1097 (Amor- 
isco et al. 2015); ram pressure stripped gas, such as the spec- 
tacular 60 kpc x 1.5 kpc Ha feature associated with the Coma 
galaxy D100 (Cramer et al. 2019); and "superthin" edge-on 
galaxies (Matthews et al. 1999). 

A gravitational lensing origin is ruled out by the identi- 
cal redshift of the galaxy that the feature points to. Tidal 
effects, ram pressure stripping, or a superthin galaxy might 
explain aspects of the main linear feature but are not con- 
sistent with the shocked gas and lack of rest-frame optical 
continuum emission on the other side of the compact galaxy. 
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Given the linearity of the entire system, the symmetry with 
respect to the nucleus, the presence of shocked gas without 
continuum emission, as well as the brightness of both the en- 
tire feature and the [O III] emission at the tip, the most viable 
explanations all involve SMBHs - either through nuclear ac- 
tivity or the local action of a set of runaway SMBHs. 


6.2. An optical jet? 


Visually, the closest analog to the linear feature is the 
famous optical jet of the z = 0.16 quasar 3C 273 (Oke & 
Schmidt 1963; Bahcall et al. 1995): its physical size is in the 
same regime (about half that of our object) and it has a simi- 
lar axis ratio and knotty appearance. However, the detection 
of bright emission lines along the feature is strong evidence 
against this interpretation. The spectra of jets are power laws, 
and there are no optical emission lines associated with optical 
jets or hot spots (Keel & Martini 1995). 

Furthermore, the 3C 273 jet and 3C 273 itself are very 
bright in the radio and X-rays, with different parts of the 
jet showing low- and high-energy emission (see Uchiyama 
et al. 2006). We inspected the VLA Sky Survey (VLASS; 
Lacy et al. 2020) as well as a 60ks deep Chandra image! of 
the field that was obtained in 2005 in the context of program 
5910 (PI Irwin). There is no evidence for a detection of the 
linear feature or the galaxy, with either the VLA or Chandra. 
We note that the z = 0.96 feature might be expected to have 
an even higher X-ray luminosity than 3C 273 if it were a jet, 
as the contribution from Compton-scattered CMB photons 
increases at higher redshifts (see Sambruna et al. 2002). 


6.3. Jet-induced star formation? 


Rather than seeing direct emission from a jet, we may be 
observing jet-induced star formation (Rees 1989; Silk 2013). 
There are two well-studied nearby examples of jets triggering 
star formation, Minkowski's object (Croft et al. 2006) and an 
area near a radio lobe of Centaurus A (Mould et al. 2000; 
Crockett et al. 2012). There are also several likely cases in 
the more distant Universe (Bicknell et al. 2000; Salomé et al. 
2015; Zovaro et al. 2019). The overall idea is that the jet 
shocks the gas, and if the gas is close to the Jeans limit sub- 
sequent cooling can lead to gravitational collapse and star 
formation (see, e.g., Fragile et al. 2017). The presence of 
both shocks and star formation along the feature is qualita- 
tively consistent with these arguments (see Rees 1989). 

The most obvious problem with this explanation is that 
there is no evidence for nuclear activity in our object from the 
BPT diagram, the VLASS, or Chandra imaging (see above). 
Itis possible, however, that the AGN turned off between trig- 
gering star formation and the epoch of observation, qualita- 
tively similar to what is seen in Hanny's Voorwerp and simi- 
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lar objects (Lintott et al. 2009; Keel et al. 2012; Smith et al. 
2022). 

A more serious issue is that the morphology of the feature 
does not match simulations or observations of jet-induced 
star formation. First, as can be seen most clearly in the top 
right panel of Fig. 1, the feature is narrowest at the tip rather 
than the base. By contrast, for a constant opening angle a jet 
linearly increases its diameter going outward from the host 
galaxy, reaching its greatest width at the furthest point (as il- 
lustrated by HST images of the M87 jet, for instance; Biretta 
et al. 1999). Second, the interaction is most effective when 
the density of the jet is lower than that of the gas, and the 
shock that is caused by the jet-cloud interaction then propa- 
gates largely perpendicular to the jet direction (e.g., Ishibashi 
& Fabian 2012; Silk 2013; Fragile et al. 2017). This leads to 
star formation in a broad cocoon rather than in the radial di- 
rection, as shown explicitly in the numerical simulations of 
Gaibler et al. (2012). It is possible for the jet to subsequently 
break out, but generically jet-cloud interactions that are able 
to trigger star formation will decollimate the jet. 

A related problem is that the observed velocity dispersion 
of the shocked gas is low. From the high resolution LRIS 
spectrum we find a velocity dispersion of < 20 km s™! in the 
main shock at the tip of the feature, which can be compared 
to ø ~ 130kms ! in the shocked gas of Centaurus A (Gra- 
ham 1998) and ø ~ 50 km s^! predicted in recent simulations 
(Mandal et al. 2021). Most fundamentally, though, the fea- 
ture is the inverse of what is expected: the strongest interac- 
tions should not be at the furthest point from the galaxy but 
close-in where the ambient gas has the highest density, and 
the feature should not become more collimated with distance 
but (much) less. 


6.4. Runaway supermassive black holes 


This brings us to our preferred explanation, the wake of a 
runaway SMBH. The central argument is the clear narrow tip 
of the linear feature, which marks both the brightest optical 
knot and the location of very bright [O 111] emission, com- 
bined with the apparent fanning out of material behind it (as 
can be seen in the top right panel of Fig. 1). As discussed be- 
low (§ 6.4.2) this scenario can accommodate the feature on 
the other side of the galaxy, as the wake of an escaped binary 
SMBH resulting from a three body interaction. The proper- 
ties of the (former) host galaxy can also be explained. Its 
compactness and irregular isophotes are evidence of the gas- 
rich recent merger that brought the black holes together, and 
the apparent absence of an AGN reflects the departure of all 
SMBHs from the nucleus. 


6.4.1. Mechanisms for producing the linear feature 


As discussed in $1 there have not been many studies of 
the interaction of a runaway SMBH with the circumgalac- 
tic gas, and there is no widely agreed-upon description of 
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Figure 7. Schematic illustration of the runaway SMBH scenario as an explanation of the key observed features. Panels 1—5 show a “classical” 
slingshot scenario (e.g., Saslaw et al. 1974). First, a merger leads to the formation of a long-lived binary SMBH (1.2). Then a third galaxy 
comes in (3), its SMBH sinks to the center of the new merger remnant, and this leads to a three-body interaction (4). One black hole (usually the 
lightest) becomes unbound from the other two and receives a large velocity kick. Conservation of linear momentum implies that the remaining 
binary gets a smaller velocity kick in the opposite direction. If the kicks are large enough all SMBHs can leave the galaxy (5). There can be 
2, 1 Gyr between the events in panels (2) and (3). Panels (4) and (5) happened ~ 40 Myr before the epoch of observation. The background 
of (6) is a frame from an Illustris TNG simulation (Pillepich et al. 2018), with lighter regions having higher gas density. This illustrates that 
there can be highly asymmetric flows in the circumgalactic medium, and we speculate that the SMBH at A is traveling through such a region 


of relatively dense and cold CGM (see text). 


what is expected to happen. Saslaw & De Young (1972) 
focus on the direct interaction between gas that is associ- 
ated with the SMBH with the ambient gas. They predict 
a strong bow shock which moves supersonically with the 
SMBH through the gas. The aftermath of the shock leads 
to a cooling cascade, ultimately leading to star formation in a 
wake behind the SMBH. de la Fuente Marcos & de la Fuente 
Marcos (2008) study the gravitational effect of the passage 
of a SMBH on the ambient gas. They find that small ve- 
locity kicks, of up to several tens of km s™!, are imparted on 
the gas, and that the subsequent new equilibrium can lead to 
gravitational collapse and star formation. There can be a de- 


lay between the passage of the SMBH and the triggering of 
star formation, depending on the impact parameter and the 
properties of the clouds. 

Both mechanisms may be important; we certainly see evi- 
dence for both star formation and shocks along the wake, in- 
cluding potentially a bow shock at or just behind the location 
of the SMBH itself, and conclude that the observations are at 
least qualitatively consistent with the models that exist. It is 
important to note that in these models the star formation does 
not take place in gas that was previously bound to the SMBH, 
but in the circumgalactic medium. The kinematics and metal- 
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licity of the gas therefore largely reflect its pre-existing state, 
perhaps slightly modified by the passage of the SMBH. 


6.4.2. Nature of the counter wake 


In this scenario there is only one explanation for the 
counter feature on the other side of the galaxy (assuming it is 
real), namely shocked gas in the wake of a second runaway 
SMBH. This is not as far-fetched as it may seem. When a 
third SMBH arrives in the vicinity of a pre-existing binary 
SMBH, a common outcome of the three body interaction is 
that one SMBH becomes unbound from the other two. The 
post-interaction binary can be the original one or contain the 
new arrival (Saslaw et al. 1974). In either case both the un- 
bound SMBH and the binary get a kick, in opposite direc- 
tions and with the velocity inversely proportional to the mass 
(Saslaw et al. 1974; Rees & Saslaw 1975). The counter fea- 
ture is then the wake of the most massive product of the three 
body interaction, namely the binary SMBH. 

The relative projected length of the  wakes is 
62kpc/36kpc-1.7:1. Here we used the location of ob- 
ject B to determine the length of the counter wake; using 
the location of the [O IIT] knot instead gives the same ratio. 
Although modified by their climb out of the potential well, 
this length ratio is likely not far from the velocity ratio of the 
black holes, at least if vgy >> vesc. Generally the least massive 
object is expected to escape (1.e., become unbound) from the 
other two in a three-body interaction, with the escape prob- 
ability occ Mg}; (Valtonen & Mikkola 1991). As the escaped 
SMBH has a lower mass than each of the two components of 
the binary, the velocity ratio between the single SMBH and 
the binary SMBH is then always > 2: 1, if linear momentum 
is conserved. A lower velocity ratio can work but only if the 
three SMBHs all have similar masses, for instance 4:4:3 for 
a:b1:b2, with b1 and b2 the two components of the binary. In 
a 4:4:3 three body interaction the probability that either one 
of the most massive objects escapes (leading to the observed 
1.7:1 ratio) is about the same as the probability that the least 
massive one escapes. 

We note that simulations indicate that complete ejections 
of all SMBHs from the halo are expected to be rare, occur- 
ring only in ^ 1 % of three-body interactions (Hoffman & 
Loeb 2007). The dynamics are complex, however, particu- 
larly when black hole spin, gravitational wave radiation, and 
gas flows into the center are taken into account (see, e.g., 
Escala et al. 2005; Iwasawa et al. 2006; Chitan et al. 2022). 
Along these lines, a modification of the simple slingshot is 
that the binary hardens due to the interaction with the third 
SMBH and merges, leading to a gravitational recoil kick. 
This could explain how the binary made it so far out of the 
galaxy, without the need for the three SMBHs to have near- 
equal masses. However, the direction and amplitude of the 
recoil depends on the mass ratios, spins, and relative orienta- 


tion of the binary at the time of the merger (e.g., Herrmann 
et al. 2007; Lousto & Zlochower 2011), and it seems unlikely 
that the two wakes would be exactly opposite to one another 
in this scenario. 

The counter wake is not only shorter than the primary wake 
in the observed u-band but also much fainter, which indicates 
that the shock has a lower velocity. The shock (and black 
hole) velocities are undetermined — although we will con- 
strain them in the next section — but as noted above, the ve- 
locity ratio between the wake and counterwake is likely 1.7. 
Assuming that the sound speed is similar on both sides of the 
galaxy, the far-UV luminosity of fast shocks is expected to 
scale with the velocity of the shock as Lyy cx Va ock (Dopita 
& Sutherland 1995). The expected ratio of the UV surface 
brightness of the two wakes is therefore 1.7? = 5, in excel- 
lent agreement with the observed ratio (also 5; see § 5). The 
post-shock pressure and temperature scale as ~ bus and 
are therefore a factor of ~ 3 lower in the counter wake. This 
may explain the lack of gravitational collapse and star forma- 
tion, although the local conditions of the CGM may also play 
a role (see $ 8). 


6.4.3. Locations of the SMBHs 


The “smoking gun" evidence for this scenario would be the 
unambiguous identification of the black holes themselves. 
The approximate expected (total) SMBH mass is Mgy ~ 
2 x 107 Mc, for a bulge mass of 7 x 10? Mo and assuming 
the relation of Schutte et al. (2019). The obvious places to 
look for them are A and B in Fig. 6. These are candidates 
for “hyper compact stellar systems" (Merritt et al. 2009), 
SMBHs enveloped in stars and gas that escaped with them. 
The expected sizes of HCSSs are far below the resolution 
limit of HST and the expected stellar masses are bounded by 
the SMBH mass, so of order 10? Mc — 107 Mo. 

Focusing first on A, the tip of the feature is compact but 
not a point source: as shown in the detail view of Fig. 5 there 
are several individual bright pixels with different colors em- 
bedded within the tip. The approximate brightness of these 
individual knots is F814W zz 29.5, after subtracting the local 
background. This corresponds to a stellar mass of 106 Mọ — 
107 Mo, in the right range for a HCSS. 

The complex tip of the feature coincides with very bright 
[O 111] emission, and an interesting question is whether this 
could be the equivalent of the narrow line region (NLR) of 
an AGN. If so, it is not composed of gas that is bound to 
the black hole, as in that case the velocity dispersion would 
be at least an order of magnitude higher. Instead, it would 
be a “traveling” NLR, with the accretion disk of the SMBH 
illuminating the neighboring circumgalactic medium as it 
moves through it. If the accretion disk produces enough 
hard UV photons to ionize the local CGM it should also 
emit X-rays. The empirical relation between [O tI] lumi- 
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nosity and X-ray luminosity of Ueda et al. (2015) implies 
Lx ~ 3 x 109 ergss^!, and with standard assumptions this 
correspond to ~ 40 counts in the existing 60 ks Chandra im- 
age. However, no object is detected, and we tentatively con- 
clude that it is unlikely that the SMBH at A is active. This is 
not definitive and further study is warranted: the Ueda et al. 
(2015) relation has significant scatter and the object is on the 
edge of the Chandra pointing, leading to a wide PSF and rel- 
atively poor point source sensitivity. 

We note that it is possible that the SMBH that is producing 
the shocks and star formation at location A is not located 
there, but is further than 62 kpc from the galaxy. In the de la 
Fuente Marcos & de la Fuente Marcos (2008) picture there 
is a delay between the gravitational impulse and the onset of 
star formation of about ~ 30 Myr. For a black hole velocity 
of ~ 10? km s^! this means that the SMBH may be several 
tens of kpc ahead of the feature. A careful inspection of the 
HST image shows no clear candidates for a HCSS beyond 
the tip. 

Turning now to object B, it is a point source at HST/ACS 
resolution that is clearly distinct from the shocked gas that 
constitutes the counter wake. However, at F814W = 25.3 (see 
$5) it is uncomfortably bright in the context of expectations 
for a HCSS. The stellar mass of B is ~ 3 x 10? Ma if the 
same M/L ratio is assumed as for the galaxy, an order of 
magnitude higher than the probable black hole mass. 

A possible explanation for the brightness of B is that it is 
a chance superposition of an unrelated object, and that the 
apparent termination of the counter wake at that location is 
coincidental. We show a detailed view of the areas around A 
and B in Fig. 8. The green bands indicate the locations of the 
[O 111] knots on each side of the galaxy, with the width of the 
band the approximate uncertainty. The [O 11] knot at the end 
of the counter wake appears to be 0"25 beyond B. Also, the 
angle between B and the galaxy is 4° offset from the angle 
between A and the galaxy. There is no obvious candidate 
HCSS at the expected location (marked by *X"), but that may 
be due to the limited depth of the 1+1 orbit ACS data. 

Finally, object C is a third candidate HCSS, but only be- 
cause of its symmetric location with respect to B. In some 
dynamical configurations it may be possible to split an equal- 
mass binary, with B and C the two components, or to have 
multiple binary black holes leading to a triple escape. These 
scenarios are extremely interesting but also extremely far- 
fetched, and without further observational evidence we con- 
sider it most likely that C is a chance alignment of an unre- 
lated object. 


7. MODELING 


Here we assume that the runaway SMBH interpretation is 
correct, and aim to interpret the details of the wake in the 
HST images in this context. In $7.1 we fit the seemingly 


Figure 8. Detailed view of the areas around A and B, in the 
summed F606W + F814W image. Green bands indicate the loca- 
tions of [O 111] knots in the LRIS spectrum. If B is a chance pro- 
jection along the line of sight, a hyper compact stellar system may 
be detectable near the cross in deeper data. In the vicinity of A, the 
complex interplay of shocks, star formation, and the SMBH itself 
could be investigated with high resolution IFU spectroscopy. 


random color variations along the wake and in $ 7.2 we link 
the line-of-sight velocity variation along the wake to spatial 
variations in the HST image. In both subsections we assume 
that the SMBH is currently located at position A and that it 
triggered star formation instantaneously as it moved through 
the circumgalactic gas. 


7.1. Stellar ages 


The color variation along the wake is shown in Fig. 9. 
The information is identical to that in Fig. 4, except we now 
show errorbars as well. Colors were measured after averag- 
ing the F606W and F814W images over 0745 (9 pixel) in the 
tangential direction and smoothing the data with a 0”15 (3 
pixel) boxcar filter in the radial direction. This is why some 
prominent but small-scale features, such as the blue pixel at 
r =42kpc, do not show up clearly in the color profile. Data 
at r > 58kpc are shown in grey as they are assumed to be 
affected by the SMBH itself (the candidate hyper compact 
stellar system “A” — see § 6.4.3). Data at r < 5 kpc are part of 
the galaxy and not of the wake. 

We fit the single burst stellar population synthesis models 
of Fig. 5 to the data. The three metallicities shown in Fig. 
5,Z=0, Z=-0.5, and Z = -1, were fit separately. Besides 
the choice of metallicity there are two free parameters: the 
overall dust content and the time since the SMBH was ejected 
Teject. The age of the stellar population 7! is converted to a 
position using 

$ 
r =62-62——. (1) 
Teject 
The best-fitting Z = —0.5 model has Ay = 1.1 and Teject = 
39 Myr, and is shown by the red curve in Fig. 9. The other 
metallicities gave similar best-fit parameters but much higher 


F606W — F814W 


Figure 9. Observed F606W — F814W color along the wake, af- 
ter smoothing with a 0715 boxcar filter. The red curve is a sim- 
ple stellar population with Z = —0.5, Ay = 1.1 mag, and age varying 
linearly with position along the wake. The best-fit time since ejec- 
tion is 39 Myr, corresponding to a projected black hole velocity of 
ven © 1600 km s™'. 


x? values. This simple model reproduces the main color vari- 
ation along the wake, with three cycles going from blue to 
red colors starting at r = 56kpc all the way to r = 15kpc. 
As noted earlier, these large and sudden color changes in the 
model curve reflect the complex evolution of red and blue 
supergiants, and are not due to a complex star formation his- 
tory. The red axis shows the corresponding age of the stellar 
population. 

The best-fitting Teject implies a projected black hole veloc- 
ity of vga © 1600kms !. This velocity is in the expected 
range for runaway SMBHs (e.g., Saslaw et al. 1974; Volon- 
teri et al. 2003; Hoffman & Loeb 2007), providing further 
evidence for this interpretation. Specifically, it is too high for 
outflows and too low for relativistic jets; besides hyperve- 
locity stars, which are thought to have a similar origin (Hills 
1988), runaway SMBHs are the only objects that are likely 
to have velocities in this range. 


7.2. Kinematics 


The black hole velocity of zz 1600 km s^! that we derive 
above is much higher than the observed line-of-sight ve- 
locities of gas along the wake, which reach a maximum of 
z 330 km s^! (see Fig. 1). The observed velocities reflect the 
kinematics of the circumgalactic medium: the passing black 


hole triggers star formation in the CGM behind it but does 
not drag the gas or the newly formed stars along with it. 

In this picture the gas and newly formed stars will con- 
tinue to move after the black hole has passed. The wake 
should therefore not be perfectly straight but be deflected, 
reflecting the local kinematics of the CGM. We show the 
F606W + F814W HST image of the wake in the middle left 
panel of Fig. 10, with the vertical axis stretched to emphasize 
deviations from linearity. The wake is indeed not perfectly 
straight, but shows several *wiggles" with an amplitude of 
~ 0.5 kpc. These deviations from a straight line are quanti- 
fied by fitting a Gaussian to the spatial profile at each position 
along the wake and recording the centroids. These are indi- 
cated with orange dots in the middle left panel and with black 
points with errorbars in the bottom right panel. 

The [O 111] A5007 velocity profile is shown in the top 
left panel, with the orange line a spline fit to the chang- 
ing velocity centroids along the wake. The velocity pro- 
file shows a pronounced change between 35 kpc and 40 kpc, 
where the line-of-sight velocity increases from ~ 150km s^! 
to zz 300 km s^! . There is a change at the same location in the 
spatial profile, suggesting that the deviations from a straight 
line are indeed correlated with the CGM motions. 

We model the connection between the line-of-sight ve- 
locities and the wiggles in the HST image in the following 
way. We assume that the black hole leaves the galaxy in a 
straight line with velocity vgy and that it triggers star forma- 
tion instantaneously at each location that it passes. The newly 
formed stars will move with a velocity 8vs4,, where vga, is the 
line-of-sight velocity measured from the [O 111] line and £ is 
a conversion factor between line-of-sight velocity and veloc- 
ity in the plane of the sky tangential to the wake. By the time 
that the SMBH reaches 62 kpc, the stars at any location along 
the wake r will have moved a distance 


62-r 


VBH 


d(r) = Bvgas(r) Q) 


that is, the velocity in the plane of the sky multiplied by the 
time that has elapsed since the passage of the black hole. 

AS Vgas is directly measured at all r, the only free parame- 
ter in Eq. 2 is A ! vgg. In practice there are several nuisance 
parameters: the model can be rotated freely with respect to 
the center of the galaxy, and there may be an offset between 
the line-of-sight velocity of the galaxy and that of the CGM 
at r - 0. We use the emcee package (Foreman-Mackey et al. 
2013) to fit for the black hole velocity and the nuisance pa- 
rameters. The number of samples is 1200 with 300 walkers; 
we verified that the fit converged. 

The best fit is shown by the red line in the bottom right 
panel and the bottom left panel of Fig. 10. The fit reproduces 
the spatial variation quite well, particularly when consider- 
ing that vg, is measured from data with 8x lower resolu- 
tion. The posterior distribution of 9 ! vgg is shown in the top 
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Figure 10. Connection between velocities along the wake and its morphology. Top left: [O 111] emission along the wake, with a fit to the 
velocity centroids in orange. Middle left: HST image of the wake, with stretched vertical axis to emphasize variations. The orange dots are 
centroids. Bottom right: Fit of a kinematic model to the HST centroids, based on the [O 111] velocity profile. This fit is also shown in the bottom 
left panel. Upper right: Distribution of posteriors for the black hole velocity vgn, modified by an unconstrained geometric parameter 6. For 
B ~ 0.3 we find that vgu is consistent with the value derived from the color variation along the wake. 


right panel. We find vgy = 05300300 km s^! The constraint 
comes directly from the amplitude of the wiggles: if the black 
hole velocity were lower by a factor two, twice as much time 
would have passed since the passage of the SMBH, and the 
wake would have drifted apart twice as much (~ 1 kpc in- 
stead of the observed ~ 0.5 kpc). 

Combining this result with that from § 7.1 we infer that the 
morphological deviations from a straight line and the colors 
of the wake can be simultaneously explained if 8 ~ 0.3, that 
is, if the gas velocities perpendicular to the wake are 30 96 of 
the line-of-sight velocities. The implied direction of motion 
is about 17? away from the line of sight (with an unknown 
component in the plane of the sky along the wake). 


8. DISCUSSION AND CONCLUSIONS 


In this paper we report the discovery of a remarkable linear 
feature that is associated with a galaxy at z = 0.96. Although 


the feature exhibits superficial similarities to other thin ob- 
jects, in particular the optical jet of 3C 273, close examina- 
tion shows that it is quite unique with no known analogs. 

We make the case that the feature is the wake of a run- 
away SMBH, relying on the small number of papers that have 
been written on this topic in the past fifty years (Saslaw & 
De Young 1972; Rees & Saslaw 1975; de la Fuente Mar- 
cos & de la Fuente Marcos 2008). This area could bene- 
fit from further theoretical work, particularly since these pa- 
pers propose a variety of formation mechanism for the wakes. 
Hydrodynamical simulations that model the shocks and also 
take gravitational effects into account might bring these ini- 
tial studies together in a self-consistent framework. 

Objects A and B are possible hyper compact stellar sys- 
tems (HCSSs; Merritt et al. 2009). Neither object is a clearcut 
case: object A is not a point source, and the actual HCSS 
would be one of several candidates within the main knot. Ob- 
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ject B is brighter than what might be expected for a HCSS 
(see Boylan-Kolchin et al. 2004; Merritt et al. 2009), and as 
we show in § 6.4.3 it may well be a chance superposition of 
an unrelated object. It could also be that Merritt et al. (2009) 
underestimate the mass that can be bound to the black hole 
(as they do not take the effects of gas or possible binarity of 
the SMBH into account), that the M/L ratio of B is much 
lower than what we estimate, or that the SMBH is more mas- 
sive than what we inferred from the galaxy mass. 

We show that the seemingly random color variation along 
the wake can be explained by a simple model of aging of the 
stars, beginning at the tip of the wake. In this interpretation 
the striking excursions in Fig. 9 are due to the varying dom- 
inance of blue and red supergiants.? The evolution of these 
stars is quite uncertain; turning the argument around, the data 
provide a validation of the qualitative behavior of the models 
from 1 to 30Myr. The implied velocity of the SMBH at A 
is vgy ~ 1600 km s ! and the velocity of the binary SMBH 
is vgy ~ 900 km s™! if the ejection was symmetric. These 
velocities are projected on the plane of the sky, and do not 
correspond to predicted line-of-sight velocities; the ratio be- 
tween the line-of-sight velocities should be ~ 1.7 but their 
absolute values are poorly constrained. 

Velocities in this range are also indicated by the straight- 
ness of the HST feature: as we show in $ 7.2 the feature is ex- 
pected to differentially disperse, and its morphology requires 
that it was created by a fast-moving object. A third piece of 
evidence for high speeds comes from the emission line ra- 
tios. As noted in $3.2 it is difficult to have [O rr1]/H ratios 
as high as ~ 10 unless there is a significant precursor compo- 
nent (photoionization ahead of the shock) and the shock has 
a velocity of at least ~ 500 kms"! (Allen et al. 2008). We 
can speculate that the precursor component may be partially 
responsible for the complexity of the tip of the feature: per- 
haps star formation is not only triggered behind the SMBH 
but also just in front of it. 

The shock velocity and luminosity provide a constraint 
on its spatial extent. From Eqs. 3.4 and 4.4 in Dopita & 
Sutherland (1996) with Lug ~ 2 x 1040 ergs s! and Vshock ~ 
1600kms'! we obtain an area of the shockfront of ~ 
0.2n^! kpc?, with n the density in cm?. For n < 0.1 (as ex- 
pected for circumgalactic gas, even with some gravitational 
compression) the shock should be resolved at HST resolu- 
tion, and possibly even from the ground. In this context it is 
interesting that there is some indication that the [O III] emis- 
sion is indeed resolved along the LRIS slit. Turning this argu- 
ment around, a high resolution image of the shock (in either 


? We note that there is no appreciable contribution from emission lines in the 
HST filters; in particular, the redshifted [O 111] doublet falls redward of the 


long wavelength cutoff of the F814W filter. 


[O ri1] or the rest-frame far-UV) could provide a joint con- 
straint on the shock velocity and the density of the gas. 

The measured line-of-sight velocities along the wake do 
not tell us much about the velocity of the SMBH and its 
accompanying shocks, but they do provide a pencil beam 
view of circumgalactic gas kinematics in a regime where 
we usually have very little information. We can compare 
the kinematics to general expectations for halo gas. The 
z = l stellar mass — halo mass relation implies a halo mass 
of ~ 3 x 10!! Mọ (Girelli et al. 2020) and a virial radius 
of ~ 80kpc (Coe 2010). Considering that the projected 
length of the wake is shorter than the physical length, the 
Tproj = 62 kpc wake likely extends all the way to the virial ra- 
dius. Using Wy = (GMyir/rvir)’> we have Vy, zz 130 kms", 
much lower than the observed peak line-of-sight velocity of 
the gas of zz 330kms !. This difference may be due to the 
passage of the SMBH itself; in the impulse approximation 
of de la Fuente Marcos & de la Fuente Marcos (2008), for 
example, the black hole imparts a velocity kick on the am- 
bient gas. An intriguing alternative explanation is that the 
trajectory of the SMBH intersected gas that is not in virial 
equilibrium but an outflow or an inflow. An example of such 
a structure is a cold stream that could be funneling gas to- 
ward the galaxy. Such streams have been seen in simulations 
(Kere$ et al. 2005; Dekel et al. 2009), although not yet ob- 
served. A cold stream could explain why the velocity disper- 
sion of the gas is so low, and perhaps also facilitated raising 
the density above the threshold needed for gravitational col- 
lapse. It might also explain why the line-of-sight velocity at 
the location of the “counter” [O 11] knot, on the other side of 
the galaxy, is much lower than the velocities along the pri- 
mary wake, and perhaps also why no star formation is taking 
place on that side. We illustrate this possibility in the right 
panel of Fig. 7. 

It is straightforward to improve upon the observations that 
are presented here. The main spectrum is a 30 min exposure 
with Keck/LRIS, and the exposure time for the near-IR spec- 
trum that was used to measure [N II]/Ha was even shorter, 
7.5 min. The extraordinary sensitivity of the red channel of 
LRIS enabled us to use the redshifted [O 111] A5007 line at 
Aobs = 9834 À for most of the analysis, despite the short ex- 
posure time. Deeper data, for instance from the JWST NIR- 
SPEC IFU, may show the expected broad, highly red- or 
blueshifted emission lines of ionized gas that is bound to the 
black holes themselves. Those data could also spatially re- 
solve flows, shocks, and star formation near A (see Fig. 8). 
The HST data is similarly shallow, at 1 orbit for each of the 
two ACS filters. Deep ultraviolet imaging with UVIS is par- 
ticularly interesting, as that could map the spatial distribution 
of shocked gas on both sides of the galaxy. A UVIS image 
would readily show whether the counter wake is real, and 
whether it points to B or is precisely opposite the main wake. 
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Finally, X-ray imaging could further constrain the physics of 
the shock and the absorbing hydrogen column (see Dopita & 
Sutherland 1996; Wilson & Raymond 1999), or even directly 
detect the accretion disk of one or more of the SMBHs. The 
currently available 60 ks Chandra image shows no hint of a 
detection but as it is very far off-axis, there is room for im- 
provement. 

Looking ahead, the morphology of the feature in the HST 
images is so striking that it should not be too difficult to find 
more examples, if they exist. Future data from the Nancy 
Grace Roman telescope can be searched with automated al- 
gorithms; this is the kind of task that machine learning al- 
gorithms can be trained to do (see, e.g., Lochner & Bassett 
2020). Although technically challenging, the most interest- 


ing wavelength to search in is probably the rest-frame far-UV, 
as it may include cases where the SMBH did not trigger star 
formation. Individual runaway SMBH systems are of great 
interest in their own right; furthermore, a census of escaped 
SMBHs can complement future gravitational wave measure- 
ments from LISA (Amaro-Seoane et al. 2017) for a complete 
description of SMBH evolution in — and out of — galaxy nu- 
clei. 
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